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ABSTRACT 

Recent infrared observations have demonstrated that photoevaporation driven 
by high-energy photons from the central star contributes to the dispersal of pro- 
toplanetary disks. Here, we show that photoevaporative winds should produce a 
detectable free-free continuum emission given the range of stellar ionizing photons 
and X-ray luminosities inferred for young sun-like stars. We point out that VLA 
observations of the nearby disk around TW Hya might have already detected 
this emission at centimeter wavelengths and calculate the wind electron density 
and mass flow rate. We also estimate the intensities of H radio recombination 
lines tracing the wind and discuss which ones could be detected with current 
instrumentation. The detection and profiles of these recombination lines would 
unambiguously prove our inference of free-free emission from photoevaporating 
disks like TW Hya. In addition, radio/millimeter data can help constraining 
wind parameters such as temperature and electron density that are fundamental 
in measuring mass flow rates. 

Subject headings: circumstellar matter — radio lines: stars — protoplanetary 
disks — stars: individual (TW Hya) 
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Introduction 



Photoevaporation driven by high-energy photons from the central star has been long rec- 
ognized by theorists as a plausible mechanis m to speed-up the cleari ng of protoplanetary disks 



even around low-mass sun-like stars (e.g., iHollenbach et al.ll2000l ). Recently, we identified 



a robust diagnostic for ph otoevaporation in the forbidden [Ne II] emission line at 12.81 fim 
( jPascucci fc Sterzild 120091 ). Relatively narrow (~20km/s) and slightly blueshifted (a few 
km/s) [Ne II] line profiles have been detected toward several evolved di sks, including those 
having a gap in their dust distribution, oft en called transitional disks ( jPascucci fc Sterzik 



2009 



Pascucci et al.ll201ll ; ISacco et al.ll2012l ) . Such profiles unambiguously point to unbound 



gas in a wind. Other wind tracers have been proposed and are being investigated, e.g. the 
[O I] line at 6300 A ( Ercolano Sz Owenl 120101 ) and the CO rovibrational band at 4.67/^m 
( jPontoppidan et al.ll201ll ). Because mass flow rates are very sensitive to the density and 
temperature of the wind, it is necessary to identify additional wind diagnostics to pin down 
these parameters. 

Here, we show that observations at millimeter and radio wavelengths can provide such 
diagnostics. The photoevaporative wind may be fully (EUV case)0 or partially (X-ray case) 
ionized and deflections of electrons in the wind by protons should result in continuum free- 
free emission (Sect. [2]). We point out that VLA observations of the nearby disk around 
TW Hya might have already detected the free-free emission from its photoevaporative wind 
(Sect. [3]). This emission should be accompanied by H radio recombination lines. We calculate 
the fluxes of a representative sample of transitions and show that they can be detected with 
current instruments (Sect. [5]). We conclude by discussing the broader implications of this 
study (Sect. E]). 



Free-Free Emission from Ionized Winds 



For protoplanetary disks, H + is the dominant ion in the surface layer ionized by stellar 
EUV and X-ray photons. In this case the thermal free-free volume em issivity can be fully 
characterized by a handful of parameters (see, e.g. iPadmanabhanl l2000h : 



6.8 x lO-^n 2 e T- l ' 2 e- hulkBT g ff [erg cm' 3 s^ 1 Hz' 1 ] 



where n e is the electron density, T is the gas temperature, and gff is the velocity- averaged 
Gaunt factor. This factor depends only on T and the frequency v of the observation and can 



^UV; 13.6 eV < hu < 100 eV 
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be computed analytically. Converting the emissivity into luminosity (L„) requires integrating 
e u dV, where V is the volume of the emitting gas. Assuming temperature is constant, L u is 
proportional to the integral of n 2 e dV =EMy, the volume emission measure. 



Assume that the disk surface absorbs a fraction / (~ 0.7 for EUV photons, see lHollenbach fc Gorti 



20091 ) of the stellar EUV photon luminosity $euv- I n steady state, / $euv is equal to the 
recombination rate of electrons and protons in the ionized gas, and the latter is propor- 
tional to EMy. Therefore, the free-free luminosity is linearly dependent on the EUV photon 
luminosity. Similarly, the absorbed X-rays create a hydrogen ionization rate balanced by 
recombinations, so that the free-free emission caused by Xray ionizations is linearly propor- 
tional to the locally incident X-ray photon luminosity. Here we use the fact that, including 
secondary ionizations, X-ray ionization rates (in s _1 ) are equal to the i ncident X-ray luminos - 
ity divided by 40 eV, the mean energy required for a single ionization ( IGlassgold et al.ll2004j ). 
Because X-rays penetrate deeper in the disk than EUV photons, the fraction / of photons 
absorbed by the d isk is lower. Detailed modeling suggests a value of 0.5 for a ~5,000K gas 
( iGorti et al.l 120111 ) . Assuming the EUV ionized layer has T=10,000K and the X-ray layer 
has T= 5,000 K, the free-free flux density at 3.5 cm can be written as: 



3.5cm 



2.9 x 10" 39 



EUV 



3.5cm 



2.4 x 10 
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51 
~d 



(2) 
(3) 



here d is the source distance in pc (where 51 pc is the distance of TW Hya, lMamajekll2005l ). 
$euv is measured in photons per second, and L x in erg/s. We note that the bulk of the 
X-ray-heated gas is typically at ~1, 000-2, 000 K. However, the free-free flux density decreases 
with gas temperature such that F 3 . 5cm (l,000 K)< 0.5 x F 3 . 5cm (5,000K). 

Direct measurements of $euv from young stars a re not available because EUV photons 
are easily absorbed by H in the interstellar medium. Alexander et al.l (120051 ) used emission 
lines at FUV wavelengths to derive order-of-magnitude estimates of $euv f° r A ve y° un g solar 
mass stars of 10 41 -10 44 s _1 . Based on this large range, we can expect a free-free 3.5 cm flux 
from the fully ionized layer ranging from 290 /iJy to 0.3 Jy at the distance of TW Hya. Stellar 
X-ray luminosities have been measured in several star-forming regions. The XMM survey 
of Taurus finds Lx values from 10 29 to 10 31 erg/s for solar-mass stars (IGiidel et al.l 120071 ) 
which convert into a free-free flux from 2.4 to 240 //Jy at 51 pc. Given that the extended 
VLA (EVLA) can reach sensitivities of a few //Jy at 3.5 cm in a few hours of integration, 
this simple calculation illustrates that free-free emission from a fully or partially ionized 
disk surface can be detected with current instrumentation out to the distance of nearby 
star-forming regions like Taurus. 
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3. The case of TW Hya 



TW Hya is a nearby, relatively old star surrounded by a transition al disk (e.g.ICa 



20021 ) . The exact structure and extension of the dust gap is debated ( lAkeson et al 



a dis cussion) as well as the mechanisms responsible for the inner clearing (e.g., 
201lh . 



vet et al 



2011 



for 



Gorti et al. 



Recently, we found that the disk of TW Hya is dispersing gas f rom its surface via 
photo evaporation driven by high-energy photons from the central star ( jPascucci fc Sterzik 
20091 ). The photoevaporative wind, as traced by the [Ne II] line at 12.81 /xm, originates 
beyond the radius that marks the transition between the thin and thic k dust disk, regardles s 
of whether that occurs at 1 or 4 AU, and extends out to about 10 AU ( IPascucci et al.ll201ll ). 



Having evidence for a fully or partially ionized disk layer and given its proximity, 
TW Hya is the ideal source to search for the type of free-free emission discussed in Sect. |5J 
Fig. [T] shows the long- wavelength portion of the s ource spectral energy distribution (SEP) 



covering from 0.87 mm out to 6 cm (data are from IWilner et al.l 12005 



Andrews et al. 2012 



and references therein). Emission from most disk grains at these long wavelengths is optically 
thin hence the flux depends on the wavelength as F v oc z/ 2+/3 ), where j3 is the wavelength 
dependence of the dust opacity. In the Log-Log plot shown in Fig. [T] the best fit to the 
millimeter data (from 0.87 to 3.4mm) is for a SED slope of 2.57±0.06, implying a /3 of 0.6. 
This is typical o f classical and more evolved disks in the Taurus-Auriga star-forming region 
( Ricci et al.ll2010l ). Our fit fully accounts for the observed 7 mm flux but it is lower than the 
3.5 cm flux by more than a factor of 2. The extra emission at 3.5 cm amounts to 140±40/xJy. 

A fully ionized wind by 4 — 6 x 10 40 EUV photons per second or a wind partially 
ionized by a star with Lx of 4 — 8 x 10 30 erg/s can alone account for this excess emission 
(see eqs. 2 and 3). The X-ray luminosity of TW Hya is measured to be ~ 1.5 x 10 30 erg/s 
(IBrickhouse et al.ll2010l ) meaning that X-rays contribute to only ~ 35/iJy of the excess flux, 
even in the optimistic assumption that all the X-ray heated gas is at 5,000 K. This small 
flux is within the la flux uncertainty, hence most of the measured excess emission at 3.5 cm 
must arise from the fully ionized EUV-layer in this disk. 

The $f,ttv derived f rom eq. 2 converts into a mass loss rate M wind of 2 — 3 x 10~ 10 M /yr 
( IHollenbach et al.lll994f ). The mass accretion rate of TW Hya is time variab le and literature 
value s span a large range, from clos e to the inferred Mm„d (5 x 10~ 10 M^/vr. iMuzerolle et al. 
2000 ) to more than 10 times higher ( Alencar fc Batalhall2002 ; Dupree et al.ll20l3 ). This range 



suggests that EUV-driven photoevaporation does not dominate yet over viscous accretion. 



We note that indirect measurements of $euv f° r TW Hya, including from the [Ne II 
range from 2 — 5 x 10 41 s -1 with a large uncertainty of about a factor of 5 (IHerczeg 



line, 



2007 
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Pascucci & Sterzikll2009h . Our estimates from the free- free continuum emission are consistent 
with these values within the reported errorbars. 

T he wind electron density n e can be also estimated from the free-free luminosity and 



EMy. iHollenbach fc Gortil ( 120091 ) found that for disks EUV-illuminated by the central star 



most of EMy comes from regions near the "gravitational radius" r g , where the hydrogen 
thermal speed is equal to the escape speed from the star gravitational field. This is in- 
deed what we see in the [Ne II] line, whose critical density is above the wind density: r g 
i s 6.2 AU for TW Hy a and most of the [Ne II] emission arises within a radius of 10 AU 
( IPascucci et al.l 120111). Additionally at the gravitational radius we have a vertical extent 
zhii ~ fg (IHollenbach fc Gortill2009[ ) . Based on these considerations, we adopt thii =10 AU 
and zhii =5AU and find n e ~ 10 5 cm~ 3 . 



Radio Continuum Emission from Other Mechanisms 



We briefly examine five additional mechanisms that could produce excess emission at 
centimeter wavelengths and discuss what observations are needed to discriminate among 
them. 

Collimated ionized outflows /jets from Class I and II sources present flat or positive spec- 
tral indexes (a > —0.1 with F v oc v a ) at ra dio wavelengths pointing to free- free emission (e.g., 



Anglada et al.lll998t iRodmann et al.ll2006l and Fig. [TJ. Hence, searches for free-free emission 
from photoevaporative winds should be carried out in evolved systems like TW Hya that 
have no jets /outflows detec ted in molecular, atomic, and/or ionic tracers (lAlencar fc Batalha 
20021 ; IPascucci et al.l 120 111 ). Most transitional disks belong to this category. We can esti- 
mate the jet contribution to the free-free emission by using the empirical relation between 
momentum rate in t he molecular outflow and the radio continuum luminosity (eq. 3 from 
Anglada et al.lll998l ). For TW Hya , even the largest measured mass accretion rate and ra- 
tio between outflow and accretion (IHartigan et al.lll995t IWhite &: Hillenbrand! 120041 ) would 
produce a 3.5 cm free-free flux of only ~30/xJy for a typical jet velocity of 200km/s, ~20% 
of the measured excess flux. Similarly, any free-free contribution from accretion shocks near 
the forming star is reduced for evolved disks. In the case of TW Hya it should amoun t to 
less than a few /xJy at 3.5 cm based on the calculations of iNeufeld fc Hollenbachl ( 119961 ) . 



Another source of radio emission in young stars is non-thermal emission originating 
in magnetic fields, also known as (gyro) synchrotron radiation. A few protostars (Class0/I) 
have been reported to have this type of emission based on spectral indexes a < —0.1 and 
about an order of magnitude variability, presumably due to magnetically induced flares (e.g., 
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Forbrich et al.ll2010l ). The easiest way to discriminate between free-free and gyro-synchrotron 
radiation is by obtaining radio observations at multiple wavelengths and computing the radio 
spectral index. In the case of TW Hya, where the only ra dio detection is at 3 .5 cm, it is time 
monitoring that can rule out gyro-synchrotron emission (IWilner et al.l 120051 ). 



Finally, both very large (mm-size) and very small (nm -size) grains can enhance the 
centimeter flux. In the case of TW Hya, IWilner et al.l (120051 ) showed that a population of 5- 
7mm sized-grains, containing 99.9% of the dust mass, can match the 3.5 cm excess emission. 
The remaining 0.1% of the dust mass is in grains of sizes between 0.005 and 1 fim in their 
model sug gesting a stric t ly bim odal dust distribution. At the other end of the grain size 
spectrum, iRafikov et al.l ( 120061 ) showed that nm-sized grains spinning at thermal rates can 
produce detectable electric dipole emission at A >0.6cm. This emission has a characteristic 
bell-like spectral shape which can be easily distinguished from the power-law spectra of 
free-free and synchrotron emission. 

As discussed in Sect. [5j free- free thermal emission from photoevaporating disks is opti- 
cally thin at millimeter and centimeter wavelengths, meaning a radio spectral index a (for 
the gas only) of -0.1. Thus, the combination of sensitive EVLA continuum observations at 
3.5 and 6 cm is the most straightforward way to constrain the radio slope of evolved disks and 
discriminate among the different mechanisms discussed here. However, the most direct way 
to test our inference of free-free emission from photoevaporating disks is via the detection of 
H recombination lines and associated blueshifts. 



5. Hydrogen Radio Recombination Lines 



In a region that is fully or partially ionized electrons will be captured by protons to 
a state n and undergo transitions to lower levels. We compute here the intensities of a 
representative sample of Una recombination lines at millimeter and radio wavelengths as a 
function of $euv an d L X - 

In our calculation, we consider spontaneous and internally stimulated emission (masing) 
but neglect externally stimulated emission. The line flux density Fi can be written as: 



Fi = B v (T)n HII 



b n rt + Tg 

n + r c 



'l_ e -(n+rch _ n _ 



(4) 



where B V (T) is the Planck function at the gas temperature T, VLhii is the solid angle sub- 
tended by the ionized wind, b n is the LTE departure coefficient for the upper state n, and 
r c , t*, and T; are the continuum a nd line optical depths (the latter corrected for non-LTE ef- 
fects) at the specific frequency V\ ( IBell & Seaquistlll978l ). For the LTE departure coefficients 
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we refer to ISalem fc Brocklehurstl (119791 ) for transitions n > 50 and to IWafmsleyl (ll990h for 



lower n, i.e. for lines that fall at millimeter wavelengths. 

Eq. H] simplifies greatly for the photoevaporative winds discussed in Sect. |2] because we 
find line and continuum optical depths that are << 1 at millimeter and radio wavelengths. 
This also means that any masing is not significant because of the relatively low electron 
densities and pathlengths. To corroborate this statement let us consider the so-called turn- 
over wavelength At, the wavelength at which r c = 1: At ~ 100/(T , ~ 1 " 35 rig22:#xr) where all 
units are in cgs except for zhii which is in pc. One sees that even a partially ionized wind at 
5,000 K becomes optically thick at wavelengths > 20 cm for plausible wind values (Sect. [2]). 
Hence, we can use the optically thin approximation to re-write eq. @]as: Fi = B v (T)Q.Hiib n Tf . 
The line optical depth can be written as r ; * = r c r*, hence F\ = F c r*b n where F c is the 
thermal continuum flux. Here, r* is the line-to-thermal continuum ratio in LTE at line 
center assuming thermal broadening of the line: 

r* = 2.33 x l^/\u- l ul 1 T- 1Ab Ei/E c (5) 

{yi in GHz, T in K) with Aui being the line widthB in kHz, and E j/E a the ratio of line t o 



continuum emission measure which we assume to be 0.9 following iBell fc Seaquistl (119781 ). 
Because F c is proportional to $euv and Lx (Sect. |5J) we obtain the following relations for 
the line flux densities assuming the Rayleigh Jeans approximation (which is valid for this 
hot gas even at millimeter wavelengths): 

F t = 2.1 x 10- 41 uihn $Euv[/iJy] (6) 

F t = 5.4 x 10- 31 v x b n Cjj L x \jiJy] (7) 

where v\ is in GHz and d is the distance in pc. 

The upper panel of Fig. [2] shows flux densities for a representative set of Hna recom- 
bination lines for a fully ionized EUV-wind (10,000 K gas) and a partially ionized X-ray 
wind (5,000 K gas). We have taken <3?euv = 10 41 s -1 and Lx = 10 30 erg/s respectively and 
scaled the fluxes to 51 pc. Line-to-free-free continuum ratios are shown in the panel below. 
They scale with the gas temperature as ~ y- 1 - 65 which explains why the cooler X-ray gas 
has higher line-to-free-free-continuum ratios. These ratios also increase as z/ 1 ' 1 x b n (with b n 
slowly decreasing) as we move to high frequencies/short wavelengths suggesting that mil- 
limetric transitions might be the easiest to detect. However, the dust thermal emission 



2 Ai>i — 2vi / cy/ 2ln(2)kT / rrih with c being the sound speed, fc the Boltzmann's constant, and the 
hydrogen mass 
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increases more steeply at shorter wavelengths as illustrated in the lower panel of Fig. [2] re- 
ducing the total line-to-continuum ratio. For the dust contribution we have taken here the 
mean 7 mm flux of transitional disks in Taurus scaled at 51 pc and then applied a power 
law of the form F v oc u 2+l3 with — 2, 0.6, to encompass the range of possible SED slopes 
(from an ISM-like to the average of Taurus disks to a SED dominated by large bodies with 
grey opacity). This plot illustrates that short cm wavelengths have actually the higher total 
line-to-continuum ratios. 

The EVLA could detect (but not spectrally resolve) some of these lines given the sen- 
sitivity of ~ 5 /zJy rms in about 3 hours with a bandwidth of ~500 MHz. In the millimeter 
regime, we will be able to detect Hna recombination lines only if the dust disk emission 
is spatially more extended than the free-free emission so that the line-to-continuum ratio 
can be increased at r g , where the wind emission peaks. Evolved disks with dust gaps and 
large $euv (or Lx) are the best candidates to detect H recombination lines at millimeter 
wavelengths. 



6. Conclusions and Perspectives 

This letter investigates the radio/millimeter properties of photoevaporative winds driven 
by high energy photons from the central star. We show that free-free continuum emission 
from fully (EUV case) or partially ionized (X-ray case) winds is directly proportional to the 
stellar ionizing flux/Xray luminosity. Given the inferred/measured &euv/Lx of young sun- 
like stars, centimeter wind emission should be detectable out to nearby star-forming regions. 
Other mechanisms producing centimeter emission can be ruled out with observations at 
multiple wavelengths measuring the radio spectral slope. However, the smoking gun for free- 
free emission from photoevaporative winds would come from the detection and profiles of H 
radio recombination lines. Our calculations suggest that a few of them should be detectable 
at radio wavelengths and might be also detectable in the millimeter with ALMA if the dust 
contribution can be spatially resolved out using high resolution. 

We point out that VLA observations might have already detected the free-free continuum 
emission from the nearby and photoevaporating disk of TW Hya. Taking TW Hya as a case 
study, we show how radio observations can be used to infer the ionizing luminosity reaching 
the disk, the wind temperature and electron density, and hence compute mass flow rates from 
star-driven photoevaporation. This latter parameter is essential to estimate disk lifetimes and 
unde rstand whether photoey aporation is primarily driven by stellar X-rays at lO~ 8 M /yr 
e.g.. lErcolano fc Owenll2010[ ) or by EUV photons at about a hundred time lower rate (e.g., 



Alexander et al.ll2006l ). Evolved disks with low mass accretion rates should be the prime 



-9 - 



targets to expand the analysis presented here and to identify what are the typical mass 
loss rates from star-driven photoevaporation. These empirically derived rates will enable 
evaluating the impact of photoevaporation on the dispersal of protoplanetary disks, as well 
as on planet formation and migration. 

IP acknowledges support from the National Science Foundation through the research 
grant AST0908479. 

Facilities: VLA, ALMA. 
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Long — wavelength SED of TW Hya 
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Fig. 1. — SED of the TW Hya disk from 0.87mm out to 6cm. Observed fluxes (empty 
circles) and 3a upper limits (downward triangle) are from Wilner et al. (2005) and Andrews 
et al. (2012). The dashed line is a linear fit to the millimeter fluxes between 0.87 and 3.4 mm. 
Note that this fit fully accounts for the 7 mm flux but under-predicts the 3.5 cm flux by a 
factor of 2.2. Dotted lines are free-free emission relations expected for optically thin (a —- 
0.1) and thick (a =0.4) gas passing through the observed minus dust emission at 3.5 cm. We 
also plot literature 3.5 cm fluxes from Class I (Anglada et al. 1998) and Class II/classical 
disks (Rodmann et al. 2006) scaled at the distance of TW Hya. 
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Fig. 2. — First panel: Line fluxes for a representative sample of Una recombination lines 
at centimeter and millimeter wavelengths. For the EUV-case (diamonds) we have assumed 
&euv — 10 s while for the X-ray case (x) we have used Lx = 10 30 erg/s. Note that line 
fluxes are directly proportional to these quantities. Second panel: Line center-to-free-free 
continuum flux ratios for the same set of transitions. Third panel: Line center-to-total 
continuum flux ratios for the same set of transitions. For the assumed continuum emission 
we refer to Sect. EJ Symbols are for /3 = 0.6, blue lines for (3 = 2, and red lines for (3 = 0. 



